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Abstract 

An updated analysis of all available neutrino oscillation evidence in Reactor 
(Kamland, first 145 days of data) and Solar experiments (SK day and night spectra, 
global rates from Homestake, SAGE and GALLEX) including the SNOCC and NC 
data is presented. In the framework of two active neutrino oscillations we deter- 
mine the allowed regions in neutrino parameter space, we obtain, from the Kamland 
spectral shape and global signal, the following antineutrino best solution parameters, 
Am 2 ^ = 7.7 x 10 _5 ey 2 , tan 2 Oki = 0.98. The overall effect of the measured observed 
ratio R ~ 0.6 is that the LMA region remains the only one which is still favored. 

Combining Kamland and Solar data and assuming CPT invariance, i.e. the same 
mass matrix for neutrino and antineutrinos, we obtain the following antineutrino 
best solution parameters (LMAI solution), Am 2 = 7.1 x W~ 5 eV 2 , tan 2 9 = 0.47. A 
second solution (LMAII) appears for values Am 2 = 1.5 x 10~ 4 eV 2 , tan 2 9 = 0.48. We 
determine additionally individual neutrino mixing parameters and their errors from 
fits to marginal likelihood distributions, the values are compatible with previous 
results. In both methods, x 2 minimization and marginal likelihood, the combined 
analysis of solar and Kamland data concludes that maximal mixing is not favored at 
the ~ 3d level at least. 
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1 Introduction 



Evidence of antineutrino disappearance in a beam of antineutrinos in the Kamland 
experiment has been recently presented [1]. The analysis of these results [1,2] in 
terms of neutrino oscillations have largely improved our knowledge of neutrino mixing 
on the LMA region. The results appear to confirm in a independent way that the 
observed deficit of solar neutrinos is indeed due to neutrino oscillations. The ability 
to measure the LMA solution, the one preferred by the solar neutrino data at present, 
"in the lab" puts Kamland in a pioneering situation: after these results there should 
remain little doubt of the physical reality of neutrino mass and oscillations. 

The publication of the SNO results [3, 4] has already made an important break- 
through towards the solution of the long standing solar neutrino [5-7] problem 
(SNP) possible. These results provide the strongest evidence so far (at least un- 
til Kamland improves its statistics) for flavor oscillation in the neutral lepton sec- 
tor. From the combined analysis of SNO and other solar evidence, one obtains, in 
the framework of two active neutrino oscillations, the following set of parameters, 
Am 2 solar = 4.5+^ x W- 5 eV 2 ,tw 2 9 so i ar = 0.40lg;Jg. We will see in this work how 
the evidence from the Kamland measurements improves or modifies these values. 

The previous generation of reactor experiments (CHOOZ [8], Palo Verde [9]), 
performed with a baseline of about 1 km. They have attained a sensitivity of 
Am 2 < 10" 3 eV 2 [8, 10] and, not finding any disappearance of the initial flux, they 
demonstrated that the atmospheric neutrino anomaly [11] is not due to muon-electron 
neutrino oscillations. The Kamland experiment is the successor of such experiments 
at a much larger scale in terms of baseline distance and total incident flux. This 
experiment relies upon a 1 kton liquid scintillator detector located at the old, en- 
larged, Kamiokande site. It searches for the oscillation of antineutrinos emitted by 
several nuclear power plants in Japan. The nearby 16 (of a total of 51) nuclear power 
stations deliver a V e flux of 1.3 x 10 6 cm _2 s _1 for neutrino energies E u > 1.8 MeV 
at the detector position. About 85% of this flux comes from reactors forming a well 
defined baseline of 139-344 km. Thus, the flight range is limited in spite of using 
several reactors, because of this fact the sensitivity of Kamland increases by nearly 
two orders of magnitude compared to previous reactor experiments. 

The aim of this work is to study the implications of the recent Kamland results 
on the determination of the neutrino oscillation parameters, to understand which 
regions of the parameter space still allowed by the solar neutrino experiments are 
favored by them. The structure of this work is the following. In section 2 we discuss 
the main features of Kamland experiment that are relevant for our analysis: The 
salient aspects of the procedure we are adopting and the results of our analysis are 
presented and discussed in section 3. Finally, in section 4 we draw our conclusions 
and discuss possible future scenarios. 
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2 The computation of the expected signals 



2.1 The Kamland signal 

Electron antineutrinos from nuclear reactors with energies above 1.8 MeV are mea- 
sured in Kamland by detecting the inverse /3-decay reaction V e + p — > n + e + . The 
time coincidence, the space correlation and the energy balance between the positron 
signal and the 2.2 MeV 7-ray produced by the capture of a already-thermalized neu- 
tron on a free proton make it possible to identify this reaction unambiguously, even 
in the presence of a rather large background. 

The two principal ingredients in the calculation of the expected signal in Kamland 
are the reactor flux and the antineutrino cross section on protons. A number of short 
baseline experiments (See Ref. [12] and references therein) have previously measured 
the energy spectrum of reactors at distances where oscillatory effects have been shown 
to be nonexistent. They have shown that the theoretical neutrino flux predictions 
are reliable within 2% [13]. The effective flux of antineutrinos released by the nuclear 
plants is a rather well understood function of the thermal power of the reactor and the 
amount of thermal power emitted during the fission of a given nucleus, which gives 
the total amount, and the isotopic composition of the reactor fuel which gives the 
spectral shape. Detailed tables for these magnitudes can be found in Ref. [12]. For a 
given isotope the energy spectrum can be parametrized by an exponential expression 
[14] where the coefficients depend on the nature of the fissionable isotope (see Ref. [12] 
for explicit values). Along the year, between periods of refueling, the total effective 
flux changes with time as the fuel is expended and the isotope relative composition 
varies. We take the average of the relative fission yields over the live time as given 
by the experiment: 235 U = 57%, 238 U = 7.8%, 239 Pu = 30%, 241 Pu = 5.7%. 

In order to obtain the expected number of events at Kamland, we sum the ex- 
pectations for all the relevant reactor sources weighting each source by its power and 
distance to the detector (table II in Ref. [12] ), assuming the same spectrum origi- 
nated from each reactor. We sum over the nearby power reactors, we neglect farther 
Japanese and Korean reactors and even farther rest-of-the-world reactors which give 
only a minor additional contribution. The average number of positrons TVj which are 
detected per visible energy bin AEi is given by the convolution of different quantities: 
P, the oscillation probability averaged over the distance and power of the different 
reactors. Expressions for the antineutrino capture cross section are taken from the 
literature [14, 15]. The matrix element for this cross section can be written in terms 
of the neutron half-life, we have used the latest published value t±/2 = 613.9 ± 0.55 
[16]. The antineutrino flux spectrum, the relative reactor-reactor power normaliza- 
tion which is included in the definition of P and the energy resolution of Kamland 
are used in addition. We use in our analysis the following expression for the energy 
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resolution in the prompt positron detection 

a{E e ) = 0.0062 + 0.065^E e . (1) 

This expression is obtained from the raw calibration data presented in Rcf. [17]. Note 
that we prefer to use this expression instead of the much less accurate one given in 
Ref.[l]. Moreover, we assume a 408 ton fiducial mass and standard nuclear plant 
power and fuel schedule, we take an averaged, time- independent, fuel composition 
equal for each detector as given above. Detection efficiency is taken close 100% and 
independent of the energy [1]. 

We will not add any background events as they (a total of 0.95 ±0.99 events which 
include random coincidences from radioactive decays, < 0.01 evts, and correlated 
background from cosmic ray muons and neutrons, < 0.5 events,) can be distinguished 
from the signal with sufficiently high efficiency or are negligible above the 2.6 MeV 
analysis threshold (table I in Ref. [1]). We also consider negligible the background 
from geological neutrinos above the 2.6 MeV analysis threshold. 

2.2 The Solar Signal 

We also need the expected signals in the different solar neutrino experiments. These 
are obtained by convoluting solar neutrino fluxes, sun and earth oscillation probabili- 
ties, neutrino cross sections and detector energy response functions. We closely follow 
the same methods already well explained in previous works [5, 18-20], we will mention 
here only a few aspects of this computation. We determine the neutrino oscillation 
probabilities using the standard methods found in literature [21], as explained in 
detail in [18] and in [5]. We use a thoroughly numerical method to calculate the 
neutrino evolution equations in the presence of matter for all the parameter space. 
For the solar neutrino case the calculation is split in three steps, corresponding to 
the neutrino propagation inside the Sun, in the vacuum (where the propagation is 
computed analytically) and in the Earth. We average over the neutrino production 
point inside the Sun and we take the electron number density n e in the Sun by 
the BPB2001 model [31]. The averaging over the annual variation of the orbit is 
also exactly performed. To take the Earth matter effects into account, we adopt a 
spherical model of the Earth density and chemical composition [22]. The joining of 
the neutrino propagation in the three different regions is performed exactly using an 
evolution operator formalism [21]. The final survival probabilities are obtained from 
the corresponding (non-pure) density matrices built from the evolution operators in 
each of these three regions. 

In this analysis in addition to night probabilities we will need the partial night 
probabilities corresponding to the 6 zenith angle bin data presented by SK [23]. 
They are obtained using appropriate weights which depend on the neutrino impact 



4 



parameter and the sagitta distance from neutrino trajectory to the Earth's center, 
for each detector's geographical location. 

3 Analysis and Results 

In order to study power of the Kamland results for resolving the neutrino oscillation 
parameter space, we have developed two kind of analysis. In the first case we deal 
with the Kamland measured global signal. In the second case we include the full 
Kamland spectrum information. We perform a complete x 2 statistical analysis before 
and after including in addition the up-date solar evidence obtaining the regions in 
the parameter space which are favored. 

3.1 Analysis of the Global rate 

The total x 2 value is given by the sum of two distinct contributions, that is the 
one coming from all the solar neutrino data and the contribution of the Kamland 
experiment: 

2 _ 2 | 2 
X — X& T Xglob,KLi 

with 

2 (R ex P - R teo (Am 2 ,e) 

Xglob,KL — I ~ 

\ o stat+sys 

The "experimental" signal ratio is R exp = 0.611 ± 0.085 ± 0.041 [1], where the first 
error is statistical and the second one systematic. 

The solar neutrino contribution can be written in the following way: 

X% = Xglob +XSK +XSNO- ( 3 ) 

The function Xgiob correspond to the total event rates measured at the Homes- 
take experiment [24] and at the gallium experiments SAGE [25,26], GNO [27] and 
G ALLEX [28] . We follow closely the definition used in previous works (see Ref . [5] for 
definitions and Table (1) in Ref. [5] ' for an explicit list of results and other references). 

The contribution to the x 2 from the SuperKamiokande data (x|k) nas Deen ob- 
tained by using double-binned data in energy and zenith angle (see table 2 in Ref. [23] 
and also Ref. [29] ) : 8 energy bins of variable width and 7 zenith angle bins which in- 
clude the day bin and 6 night ones. The definition is given by: 

xIk = (^-^^(aL + aL)- 1 ^-^). (4) 

The theoretical and experimental R quantities are this time matrices of dimension 
8x7. The covariance quantity a is a 4-rank tensor constructed in terms of statistic 
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errors, energy and zenith angle bin-correlated and uncorrelated uncertainties. The 
data and errors for individual energy bins for SK spectrum has been obtained from 
Ref. [23]. 

The contribution of SNO to the \ 2 1S given by 

X 2 pec-SNO = Y.^-^yi^t+^ysty'iR^-R^), (5) 

d,n 

where the day and night R vectors of dimension 17 are made up by the values of the 
total (NC+CC+ES) SNO signal for the different bins of the spectrum. The statistical 
contribution to the covariance matrix, cr s t a t, is obtained directly from the SNO data. 
The part of the matrix related to the systematical errors has been computed by us 
studying the influence on the response function of the different sources of correlated 
and uncorrelated errors reported by SNO collaboration (see table II of Ref. [4]), we 
assume full correlation or full anticorrelation according to each source. 

To test a particular oscillation hypothesis against the parameters of the best fit 
and obtain allowed regions in parameter space we perform a minimization of the two 
dimensional function x 2 (Am 2 , tan 2 9). A given point in the oscillation parameter 
space is allowed if the globally subtracted quantity fulfills the condition A% 2 = 
X 2 (Am, 2 ,#) - Xmi n < Xn{CL). Where x 2 =2( 90 %, 95%, ...) are the quantiles for two 
degrees of freedom. 

In Figs.l we graphically show the results of this analysis. The favored region 
by the Kamland global rate alone is presented in fig. 1 (Left). Fig. 1 (right) includes 
the full Kamland and Solar contributions. Kamland global rate alone restricts the 
value of tan 2 to be in the range tan 2 9 > 4 x 10 -1 — 1 (symmetric with respect 
the value 6 = ir /4) except for some small region corresponding to tan 2 9 ~ 1 and 
Am 2 ~ 2 — 3 x 10~ 5 eV 2 . The overall effect of the measured observed ratio R ~ 0.6% 
is that the LMA region remains the only one which is still favored. Only if Kamland 
should had measured a ratio ~ 0.1 — 0.3, more than three sigma away from the actual 
measurement, other regions as the LOW solar solution would have some statistical 
chance to survive. In Table (La) we present the best fit parameters or local minima 
obtained from the minimization of the \ 2 function. Also shown are the values of 
Xmin- According to these results, regardless of the individual Kamland results, a 
maximal mixing solution is excluded when solar data is included in the analysis. 

3.2 Analysis including the full Kamland signal 

Here we use the binned Kamland signal (See table 2 extracted from Fig. 5 in Ref.[l]) 
combined with the evidence of the up-date solar experiments (CL,GA,SK and SNO). 
The total \ 2 value is given now by the sum of two distinct contributions, that is the 
one coming from all the solar neutrino data and the contribution of the Kamland 
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experiment which includes both, shape and total signal information: 

2 2,2 | 2 f R \ 

X — Xq t Xspec,KL "+~ Xglob,KL- \°) 

The contributions of the solar neutrino experiments and global Kamland signal 
are described in detail in the previous section. The contribution of the Kamland 
spectrum is now as follows: 

X s 2 pec,KL = (aR^-R^rial^ + al^iaR^-R^) (7) 

The total error matrix a is computed as a sum of assumed systematic deviations, 
Csys/S ~ 6.5%, mainly coming from flux uncertainty (3%), energy calibration and 
threshold (see table II of Ref. [1] for a total systematic error ~ 6.4%), see also Ref.[12, 
17,30]) and statistical errors. The parameter a is a free normalization parameter. 
The effect of systematic sources on individual bin deviations has been computed by 
us studying the influence on the response function, furtherly we have assumed full 
correlation among bins. 

Note that, as an alternative to the previous formula 7 based on ratios and using 
and standard \ 2 expression, we could have directly used the absolute number of events 
in each bin and applied poissonian statistics formulas. The Kamland collaboration 
in its published analysis uses a mixed strategy and divides the x 2 expression in 
two parts: it adds a standard \ 2 contribution corresponding to the global rate to 
a likelihood contribution which models the statistical significance of the spectrum 
shape. In this work, where we are interested not only on the Kamland data itself but 
in combining it with the solar data, we prefer to continue using the familiar expression 
7 for various reasons: first, for simplicity, clarity and comparability with previous and 
future results, second, in the major part of the bins, certainly in the most significant 
ones, the signal is high enough, and more importantly the errors are small enough, 
to use the Gaussian approximation. For the large energy bins the signal is strictly 
zero: its contribution is anyway effectively zero in both approaches. Finally, not less 
importantly, the possibility of introducing in our approach in a straightforward way 
the effect of correlated systematic deviations among bins. 

The R are length 13 vectors containing the binned spectrum (0.425 MeV bins 
ranging from 2.6 to 8.125 MeV) normalized to the non-oscillation expectations. The- 
oretical vectors are a function of the oscillation parameters: R th = R th (Am 2 ,9). 
The experimental vectors _R cxp contain the Kamland binned signal. We generate 
acceptance contours (at 90,95 and 99 % CL) in the (Am 2 , tan 2 6) plane in a similar 
manner as explained in the previous section including now a minimization respec the 
parameter a for any of the other oscillation parameters. For the sake of compar- 
ison we have also obtained exclusion regions derived from the consideration of the 
Kamland contribution alone (x 2 KL = X%ec + xli)- 
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In Figs. (2) we graphically show the results of this analysis. The first case, study 
of the Kamland data alone (Xkl) is represented by the Fig.(2)( left). The allowed 
regions in parameter space corresponding to each particular point is formed by a 
number of regions symmetric with respect the line tan 9 = 1. The general effect of 
the inclusion of the solar evidence, showed in Fig.2 (right), in the x 2 is the breaking 
of the symmetry in tan 2 9, as expected, and the strong contraction of the allowed 
area to well defined regions. The best fit parameters obtained from the minimization 
of the x 2 function are presented in Table (La) where two well separated solutions 
appear LMAI,LMAII. In both cases the mixing angle is around tan 2 9 ~ 0.45 — 0.48. 

Again, the introduction of the solar data strongly diminishes the favored value for 
the mixing angle. The final value is more near to those values favored by the solar 
data alone than to the Kamland ones. This effect could be simply due to the present 
low Kamland statistics or, more worrying, to some statistical artifact derived from 
the complexity of the analysis and of the heterogeneity of binned data involved. 

In this respect, we perform additionally a second kind of analysis in order to 
obtain concrete values for the individual oscillation parameters and estimates for 
their uncertainties. We study the marginalized parameter constraints where the x 2 
quantity is converted into likelihood using the expression C = e - ^ 2- ^™™^ 2 . This 
normalized marginal likelihood, obtained from the integration of C for each of the 
variables, is plotted in Figs. (3) for each of the oscillation parameters Am 2 and 
tan 2 9. For tan 2 9 we observe that the likelihood function is strongly peaked in a 
region tan 2 9 ~ 0.4 — 0.6. The situation for Am 2 is similar except for the existence of 
additional, somehow narrower, secondary peaks. Concrete values for the parameters 
are extracted by fitting one- or two-sided Gaussian distributions to any of the peaks 
(fits not showed in the plots). In both cases, for angle and the mass difference 
distributions the goodness of fit of the Gaussian fit to each individual peak is excellent 
(g.o.f > 99.8%) thus justifying the consistency of the procedure. The values for the 
parameters obtained in this way appear in Table (Lb). They are fully consistent and 
very similar to the values obtained from simple x 2 minimization. In particular, the 
maximal mixing solution is again excluded at the ~ 3u level. 

Although both are mutually compatible, the slight difference of the value obtained 
for the mixing angle is well explained by the shape of the allowed regions in Fig 2 
(right): the right elongation of these makes the value of the integral which defines 
the marginal distribution for tan 2 9 to be shifted. Additional variability can be eas- 
ily introduced if would have used different prior information or different parameter 
definition (i.e. sin 29 instead of tan 2 9). This an example of how the details of statisti- 
cal analysis can significantly modify the values of the physical parameters extracted 
from data and how important is that different methods are explored, specially in 
the present context of neutrino physics where heterogeneous data coming from very 
different sources is jointly used. 
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4 Summary and Conclusions 



We have analyzed the present experimental situation of our knowledge of the neutrino 
mixing parameters in the region of the parameter space that is relevant for solar 
neutrinos and we have studied in detail how this knowledge has improved with the 
recent results presented by the reactor experiment Kamland. We show, how, in 
general, the regions selected by Kamland alone, all symmetric with respect to tan 2 9 = 
1, have a large spreadth in the mixing angle. The experiment has, however, a much 
higher sensitivity to the mass difference parameter. 

Kamland global rate alone restricts the value of tan 2 9 to be in the range tan 2 9 > 
4 x 10 -1 — 1 (symmetric with respect the value 9 = vr/4) except for some small region 
corresponding to tan 2 9 ~ 1 and Am 2 ~ 2 — 3 x 1CP 5 eV 2 . The overall effect of the 
measured observed ratio R ~ 0.6% is that the LMA region remains the only one 
which is still favored. Only if Kamland should had measured a ratio ~ 0.1 — 0.3, 
more than three sigma away from the actual measurement, other regions as the LOW 
solar solution would have some statistical chance to survive. 

From the analysis of the Kamland evidence alone, the allowed regions in parameter 
space corresponding to each particular point is formed by a number regions symmetric 
with respect the line tan 9 = 1. The general effect of the inclusion of the solar evidence 
in the analysis is the breaking of the symmetry in tan 2 9, the strong reduction of the 
area of the allowed regions and, more importantly, the shifting of the value of the 
best fit mixing angle from being compatible with maximal mixing to more than 3a 
incompatibility. 

In addition to parameter extraction from \ 2 minimization, we obtain concrete 
values for the individual oscillation parameters and estimates for their uncertainties 
from the marginalized likelihood distributions. For tan 2 9 we observe that the likeli- 
hood function is concentrated in a region slightly higher that that one obtained from 
minimization. According to this analysis, maximal mixing is again excluded at the 3a 
level. The marginal distribution for Am 2 shows clearly the existence of an additional, 
somehow narrower and much less significant, secondary peaks. Kamland, after only 
150 days of data taking, together with the rest of solar experiments is already able to 
resolve the neutrino mass difference with a high precision of SAm1 l+solar < ±0.8 to be 
compared with the solar only case where the precision obtained was 5Am 2 olar < ±1.5. 
The situation for the other oscillation parameter is worse, from the present data we 
obtain 5 tan 2 9ki+ so iar < ±0.15. This precision is comparable, in fact even worse, 
than the precision attainable from the solar data alone 5 tan 2 9 so i ar < ±0.10 [5]. The 
perspectives of Kamland, after 1-3 years of data taking, together or not with solar 
evidence, for a much better determination of the mixing angle are not very optimistic 
(see the results in Ref.[15]). The ability of experiments as BOREXINO to improve 
the determination of the mixing angle with a smaller error in the near future remains 
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an open question (see Ref.[32]). 

In summary, in the framework of two active neutrino oscillations we obtain, com- 
bining Kamland and Solar data and assuming CPT invariance, i.e. the same mass 
matrix for neutrino and antineutrinos, the following values for neutrino mixing pa- 
rameters Am^ = 7.1 x W~ 5 eV 2 , tan 2 9ki = 0.47 from \ 2 minimization (LMAI solu- 
tion). 
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Am 2 (eV 2 ) tan 2 x 2 



Table La (Minimization \ 2 )'- 

KL (Sp) 7.7 x 1(T 5 0.98 1.93 

KL (Gl)+Solar 5.8 x lO" 5 0.47 45.8 

KL (Sp+Gl)+Solar, LMAI 7.1 x 10~ 5 0.47 49.1 

KL (Sp+Gl)+Solar, LMAII 1.5 x 10~ 4 0.48 49.8 

Table Lb ( From Fit, Fig.3, ±la): 

KL (SP+Gl)+Solar, LMAI 8.0±g;l x lO" 5 

KL (SP+Gl)+Solar, LMAII 1.7±g;l x 10~ 4 

KL (SP+Gl)+Solar 0.55+°;}° 

Table 1: Mixing parameters: From \ 2 minimization (Tables La) and from fit to the peak of 
marginal likelihood distributions (Table Lb). 



Bin (MeV) 


R — S exp / Smc 




2.600-3.025 


0.435 


±0.160 


3.025-3.450 


0.689 


±0.215 


3.450-3.875 


0.666 


±0.225 


3.875-4.300 


0.719 


±0.250 


4.300-4.725 


0.885 


±0.310 


4.725-5.150 


0.550 


±0.305 


5.150-5.575 


1.000 


±0.460 


5.575-6.000 


0.598 


±0.500 


6.000-6.425 


0.000 


±0.365 


6.425-6.850 


0.000 


±0.630 


6.850-7.275 


0.000 


±2.500 


7.275-7.700 


0.000 


±2.500 


7.700-8.125 


0.000 


±2.500 



Table 2: Summary of Kamland signal ratios as used in the present work. The table has 
been obtained from the information contained in Fig.5 in Ref.[l]. 
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Figure 1: Exclusion plots including Kamland global rates. The colored areas are allowed at 90, 
95, 99 and 99.7% CL relative to the absolute minimum. The region above the upper thick line is 
excluded by the previous CHOOZ experiment [10]. 




Figure 2: Allowed areas in the two neutrino parameter space after 150 days of data taking in 
Kamland. The colored lines separate allowed regions at 90, 95, 99 and 99.7% CL relative to 
the absolute minimum. (Left) Results with the Kamland signal alone (shape and total signal 
information, only 90,95% CI lines are shown for clarity). (Right) Kamland spectrum plus solar 
(CL,GA,SK,SNO) evidence. 



15 







0.2 



0.5 1 
tan 2 6 




10~ 5 10~ 4 
Am 2 /eV 2 

Figure 3: Marginalized likelihood distributions for each of the oscillation parameters Am 2 (right), 
tan 2 6 (left) corresponding to the Kamland Spectrum plus solar evidence (Fig.2(Right)). The 
curves are in arbitrary units with normalization to the maximum height. Values for the peak 
position are obtained by fitting two-sided Gaussian distributions (not showed in the plot). Dashed 
lines delimit ±lcr error regions around the maximum. See Table La for values of the position and 
widths of the peaks. 
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